Breit Interaction and Parity Non-conservation in Many-Electron Atoms 
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We present accurate ab initio non-perturbative calculations of the Breit correction to the parity 
non- conserving (PNC) amplitudes of the 6s — 7s and 6s — 5(^/2 transitions in Cs, 7s — 8s and 7s — &d 3 /2 
transitions in Fr, 6s — 5d 3 / 2 transition in Ba + , 7s — 6d 3 / 2 transition in Ra + , and 6pi/ 2 — 6p 3 / 2 
transition in Tl. The results for the 6s — 7s transition in Cs and 7s — 8s transition in Fr are in good 
agreement with other calculations while calculations for other atoms/transitions are presented for 
the first time. We demonstrate that higher-orders many-body corrections to the Breit interaction 
are especially important for the s — d PNC amplitudes. We confirm good agreement of the PNC 
measurements for cesium and thallium with the standard model . 

PACS numbers: PACS: 32.80.Ys,31.30.Jv,11.30.Er,12.15.Ji 



I. INTRODUCTION 

Study of the parity non-conservation (PNC) in atoms 
have reached the accuracy where small corrections like 
Breit interaction and radiative corrections play signifi- 
cant role (see, e.q. review pj). For example, experimen- 
tal accuracy for the PNC in cesium is 0.35% Q|. Inter- 
retation of this result using earlier atomic calculations 
~| lead to apparent disagreement with the standard 
model GaL This disagreement had been resolved when 
Breit 0, 0, H| and radiative 0, 0, |H| corrections were 
included into analysis. 

In present paper we revisit Breit corrections to the 
PNC amplitudes in many-electron atoms. Previous cal- 
culations were focused on the 6s — 7s PNC amplitude in 
Cs [111 SEI, while some also considered very similar 
(in terms of electron structure) 7s — 8s PNC amplitude 
in Fr 0, |2(j. The results are considered to be in 
reasonable agreement with each other in spite of some 
difference in the final values given by different calcula- 
tions (up to ~ 30%). This is because Breit correction to 
the PNC amplitude in cesium was small (~ —0.6%) and 
relatively rough estimations were sufficient for the sake 
of the interpretation of the experimental measurements. 

However, if we need to calculate Breit corrections to 
PNC in atoms heavier than cesium or, in general, if we 
want to have a reliable method to study the effect of Breit 
interaction in heavy many-electron atoms the important 
question to ask is what kind of approach should be used 
to do this? The answer is evident from the analysis of 
the details of the calculations of Breit contributions by 
means of many-body perturbation theory (MBPT). The 
analysis reveals very poor convergence of the MBPT with 
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respect to the number of residual Coulomb interactions 
included into the higher-order terms. This means that 
some all-order technique is needed to sum up important 
chains of higher-order terms to all orders. The natural 
choice is to present inter-electron interaction as a sum of 
Coulomb and Breit terms and to treat both terms the 
same way. 

We have developed such approach in our previous 
work and applied it to the analysis of the PNC in 
cesium. Breit and Coulomb interaction were equally 
treated on the self-consistent field level (Hartree-Fock 
and random-phase approximation calculations). How- 
ever, inclusion of the Breit interaction into calculation 
of the correlations was not comprehensive. In that work 
we included correlations by calculating the correlation 
potential E for valence electrons, using it to calculate 
the so-called Brueckner orbitals for the states of valence 
electrons, and then replacing Hartree-Fock orbitals by 
Brueckner orbitals in the EIpmc amplitude. Brueckner 
orbitals were calculated with Breit interaction included 
into self-consistent Hartree-Fock potential while correla- 
tion potential £ was calculated without inclusion of the 
Breit interaction. 

In present work we do one more step by including 
Breit interaction into second-oder correlation potential 
£ thus achieving equal and comprehensive treatment of 
Coulomb and Breit interactions on levels up to second or- 
der of the many-body perturbation theory. We apply this 
approach to the calculation of the Breit corrections to the 
PNC amplitudes in neutral cesium, thallium and fran- 
cium as well as barium and radium positive ions. This 
is the full list of atoms and ions with one external elec- 
tron above closed shells for which PNC measurements 
were done or planned. Our calculations represent the 
most comprehensive treatment of the Breit contribution 
to the PNC in many-electron atoms. The results for Tl, 
Ba + , Ra + as well as for the s — d transitions for Cs and 
Fr are presented for the first time. 

The method developed in present work is important 
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not only for the PNC in many- electron atoms. It can be 
equally useful for calculation of Breit interaction con- 
tribution to any observable values enhanced on short 
distances. These include magnetic dipole and electric 
quadrupole hyperfine structure, isotope shift, T-odd ef- 
fects, etc. 



II. THEORY 

We use the following form of the Breit operator [l^ 
(atomic units are used throughout the paper) 



jjB _ Q i' a 2 + ( Q i ' n)(«2 ■ n) 
2r 



(1) 



Here r = hr, r is distance between electrons and cti is 
a-matrix of the corresponding electron. This is a low fre- 
quency limit of the relativistic correction to the Coulomb 
interaction between electrons. It contains magnetic inter- 
action and retardation. 

Similar to Coulomb interaction, Breit interaction cre- 
ates self-consistent Breit potential 

V B 4>{r) = ]T( J^ n H B i> n d 3 r'ij(r) - 

^ n H B ^d 3 r' > n (r)). (2) 



Here summation goes over all core states. In the case 
of closed-shell atoms (like alkali atoms in the V N_1 ap- 
proximation) direct term in Breit potential vanishes (first 
term in J5J)) and only second, exchange, term remains. 
Hartree-Fock Hamiltonian has the form 

Ze 2 

H Q = cap + (/3 - l)mc 2 + V, (3) 



where V is the self-consistent potential created by elec- 
trons from the core. To achieve comprehensive treatment 
of the Breit interaction potential V should be presented 
as a sum of Coulomb and Breit terms 



V = V 



V 1 



(4) 



Note that we write for simplicity nuclear potential in © 
as for point-like nucleus. However, our actual potential 
takes into account finite nuclear size. 

To do PNC calculations one needs to take into account 
dipole interaction of an atom with external photon as 
well as weak interaction of atomic electrons with the nu- 
cleus. We use the time-dependent Hartree-Fock method 
(TDHF) to do this 16}. This method is equivalent to 
well-known random-phase approximation (RPA) . 

For every electron in the atom, single electron wave 
function is presented in the form 



tpn = Ipn + Hn + X a e~ iut 

+ 5X n e~ wt + SY n e luJt , 



Y n e l 



where index n numerates single-electron states, i/> ra is un- 
perturbed wave function for the state n which is an eigen- 
state of the Hartree-Fock Hamiltonian ® , Sip n is correc- 
tion to ip n due to weak interaction with the nucleus, X n 
and Y n are corrections due to electric field of external 
photon of frequency u>, 5X n and SY n are corrections due 
to simultaneous action of the weak interaction and dipole 
interaction with external photon. 

Corrections 8ip n , X n , Y n , SX n , SY n to all atomic states 
found by self-consistent iteration of the TDHF equations: 



(H - e n )Stjj n 
(Ho — e n — w)X n 
(Ho - e„ + cj)F„ 
(H - e„ - Lu)8X n 



(Ho 



to)SY n 



-(Hei + SV E1 )ip n , 

-6V E lSlp n 

SV W X n - SVsiWlpn, 

-5V E1 5ip n 

SV\vY n - SV Elw ip n . 



(6) 
(7) 

(8) 



Here 8V is correction to the self-consistent Hartree-Fock 
potential V due to external field (weak (W) , dipole elec- 
tric (El) or both (EIW)). Equations P7IS1) are first 
solved self-consistently for states in the core. Then cor- 
rections to valence states are calculated in the field of 
frozen core. 

PNC amplitude for the transition between valence 
states a and b in the TDHF approximation is 



jpTDH F 
^PNC 



(i> b \H E1 + 5V E1 \5ip a ) (9) 

{lpb\H W + SVw\X a ) + (lpb\SVElw\tpa). 



Breit interaction is included by using expression Q for 
the self-consistent potential on all stages of the calcu- 
lations. In other words, Hartree-Fock potential in the 
left-hand side of equations I|tjl7l8[l is the sum of Coulomb 
and Breit terms. Correction to the potential SV in the 
right-hand side of equations (|6I7I8I9H is also the sum of 
Coulomb and Breit terms: 



6V = 6V° +SV B . 

Expression does not take into account correlations. 
To include correlations we introduce correlation potential 
S. It is defined in such a way that its average value over 
particular valence state a is the correlation correction to 
the energy of this state (see, e.g. for details): 



(a\t\a). 



(10) 



(5) 



£ is a non-local operator which depends on energy and 
parameter n which defines the angular part of wave func- 
tion ip a : n a = ( — l) ia+: ' a+1/ ' 2 (j a + l/2), where l a is angular 
momentum and j a is total momentum of state a. 

In linear in S approximation, correlation correction to 
the PNC amplitude is 

AE PNC = (i/j b \£ Kb (e b )\5X a ) + {5ip b \±- Kb (e b )\X a ) 

+ (<Sy 6 |i5 Ko (e a )|^«> + (y 6 |S_ Ka (e a )|^a). (if) 
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FIG. 1: Second-order correlation potential B(ri,r2). Left arrows denote core states, right arrows denote states above the core, 
waved line is the sum of Coulomb and Breit interactions. 



In the most precise calculations of the PNC in ce- 
sium 0, 0] , operator E was calculated in all orders in 
Coulomb interaction. For the purpose of the present work 
we can restrict ourself with the lowest, second-order ex- 
pression for E. Corresponding Brueckner-Goldstone di- 
agrams are presented on Fig. ^ Breit interaction is in- 
cluded in calculation of E in two ways. First, a complete 
set of single-electron states used to calculate E was ob- 
tained with the use of the Hartree-Fock Hamiltonian © 
with Breit potential included as in eq. (@J. Second, inter- 
electron interaction in E (wave lines in all diagrams on 
Fig-© is the sum of Coulomb (e 2 /r) and Breit £Q) terms. 
The total PNC amplitude in linear in E approximation 



Epnc 



Epnc F + &-Epnc- 



(12) 



We also use a more accurate way of calculating correla- 
tion corrections which includes correlation potential E to 
all orders. We do this with the help of so-called Brueck- 
ner orbitals 16|. Brueckner orbitals are found by includ- 
ing E into Hartree-Fock equations for valence states. For 
example, for the states a and b Brueckner orbitals tp B r 
and ipa T and corresponding energies e a and e& are found 
by solving 



(#o + E 
E 



Br 



(H + E Kb - e b )^ Br 



(13) 



Corrections Sip Br , X Br and Y Br to Brueckner orbital 
ip Br induced by weak and electromagnetic interactions 
are found by solving 



(Ho + E_ 
(H + E_ K , - < 
(Ho 



-Kb 
-Kb 



-")x* r - 



£_ K .„-e a +Lo)Y a Br 



-(H w + 6V w )^ Br , 
-(H E1 + 5V E1 )iP* r , (14) 
-(H E1 +SV E1 )^ Br . 



Then, the all-order in E PNC amplitude is given by the 
expression similar to but with Hartree-Fock orbitals 
replaces by Brueckner orbitals: 



E 



Br 
PNC 



{^ r \H E1 +SV E1 \6^ r ) 



(15) 



+ (^ r \Hw + 6V w \X» r ) + W\8V E iwW) 



Equations H12|) and (|15[) differ by non-linear in E terms 
(E 2 , E 3 , etc.). 



Finally we note that non-perturbative treatment 
of Breit interaction leads to inclusion of higher- 
order in Breit operator terms, terms proportional to 
(H B ) 2 ,(H B ) 3 , etc. with H B given by (JTJ). Inclusion 
of these terms cannot be justified since expression Q 
is already approximate. Non-linear terms can be eas- 
ily eliminated by a rescaling procedure. In this proce- 
dure inter-electron interaction is considered as a sum 
U(r) = e 2 /r + XH B in which Breit interaction H B is 
multiplied by a scaling parameter A. By running pro- 
grams for different values of A a range of values is found 
for which final answer is linear function of A. Then the 
answer is interpolated to A = 1. It turns out, however, 
that as a rule A — 1 is already in linear regime. 



III. RESULTS AND DISCUSSION 

Table H] presents results of calculations of the Breit cor- 
rection to the PNC amplitudes in Cs, Ba, Tl, Fr and Ra. 
This is the full list of atoms with one external electron 
above closed shells for which PNC measurement has been 
carried out or planned. Comparison with other calcula- 
tions for the 6s — 7s transition in Cs and 7s — 8s tran- 
sition in Fr are presented in Table IVII The results for 
other atoms and for the s — d transitions in Cs and Fr 
are presented for the first time. 

Note that for all PNC amplitudes presented in Table |U 
the effect of correlations on Breit correction is small in 
E^ 2 ) approximation. It becomes significant for some s-d 
amplitudes only in Brueckner approximation. This is be- 
cause the correlation corrections in d-wave are large and 
significantly change density of the valence electron inside 
electron core. Therefore, we should perform calculations 
using Bruckner orbitals which include E to all orders. 
The linear in E expressions (lllf) . <|12|) do not provide sat- 
isfactory accuracy. 

As one can see from Table|I]the behavior of s—s and s— 
d PNC amplitudes is very different. Let us discuss them 
separately. It is convenient to use the sum-over-states 
expression for the El PNC amplitude for the discussion: 



El 



PNC 



(ab) 



E 

n 

E 



(b\D\n)(n\W\a) 

E a — E n 

(b\W\n)(n\D\a) 
Eb — E n 



(16) 
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TABLE I: Contribution of Breit interaction (A) into Epnc transition amplitudes in Cs, Ba + , Tl, Fr and Ra + (10 11 iea,BQw /N) 



Atom 



Transition 



Approximation Coulomb only Coulomb + Breit 



A/E PNC (%) 



Cs 



55 



Sl/2 



■78! 



/2 



TDHF" 

£(2)6 
Br c 



0.8903 
0.9316 
0.9001 



0.8849 
0.9258 
0.8945 



-0.0054 
-0.0058 
-0.0056 



-0.61 
-0.62 
-0.62 



6si/ 2 - 5d 3/2 



TDHF" 

£(2)6 

Br c 



3.1177 
3.5371 
3.7277 



3.1012 
3.5195 
3.7165 



-0.0165 
-0.0176 
-0.0112 



-0.53 
-0.50 
-0.30 



Ba + 56 6s 1 / 2 — 5d 3 / 2 



TDHF" 

£(2)6 

Br c 



2.1678 
2.2091 
2.1826 



2.1513 
2.1915 
2.1653 



-0.0165 
-0.0176 
-0.0173 



-0.76 
-0.80 
-0.79 



Tl 81 6p 1/a - 6p 3 /2 



TDHF a 

£(2)6 

Br c 



30.384 
23.157 
24.270 



30.105 
22.941 
24.046 



-0.279 
-0.216 
-0.224 



-0.92 
-0.93 
-0.92 



Fr 87 7s 1/2 - 8s 1/2 



TDHF a 

£(2)6 

Br c 



16.043 
15.884 
15.229 



15.891 
15.729 
15.079 



-0.152 
-0.155 
-0.150 



-0.95 
-0.98 
-0.98 



7s 1/2 - 6d 3/2 



TDHF a 

£(2)6 

Br c 



52.118 
55.568 
56.437 



51.678 
55.103 
56.090 



-0.440 
-0.465 
-0.347 



-0.84 
-0.84 
-0.61 



Ra^ 



7s!/ 2 - 6d 3 / 2 



TDHF a 

£(2)6 

Br c 



44.038 
42.817 
42.511 



43.484 
42.269 
41.971 



-0.554 
-0.548 
-0.540 



-1.26 
-1.29 
-1.27 



a TDHF: time-dependent Hartree-Fock, Eq.£jJ 

^E^ 2 ): second-order correlations are included, Ea.ll2t 

c Br: Brueckner orbitals, Eq. 1151 



where D is electric-dipole operator and W is weak in- 
teraction operator. In the RPA approximation, D = 
Hei+SVei and W = Hw+5Vw- Equation (|Fo)l is an ac- 
curate expression for the El PNC amplitude if all states 
| a), |n) are exact many-electron states of the whole 
atom. However, for our discussion we will consider only 
states which differ from the ground state by excitation 
of the single valence electron. This is reasonably good 
approximation for alkali atoms and it is sufficient for the 
discussion. In this case matrix elements in (|16|l are re- 
duced to single-electron matrix elements and we will refer 
to this expression as sum over single-electron states. 

As one can see from Table H] Breit contribution to the 
correlation correction to the s — s PNC amplitudes is 
small. This can be explained the following way. Breit and 
weak interactions are both short-range operators. They 
have significant value in the vicinity of nucleus where all 
single-electron wave functions are similar and differ by 
normalization only. Indeed, Dirac equation for all low 
states is the same on short distances since energy can 
be neglected compared to the nuclear potential. On the 
other hand, higher states with large energies do not con- 
tribute to the sum Ijltjll due to large energy denominator. 
Therefore, relative value of Breit correction to all weak 



matrix elements in <|10l) is the same and can be presented 
as a common factor outside of the summation. The effect 
of correlations on short distances can also be reduced to 
change of normalization of the wave function. Therefore, 
relative Breit correction should remain the same. 

This way of argument was used in calculatingradiative 
corrections to the PNC amplitude in cesium j^, IToL [Til 
|P2| . Moreover, it was assumed in ^(| that to find radia- 
tive correction it is sufficient to consider only corrections 
to the weak matrix elements. Radiative corrections to El 
transition amplitudes and energies were not considered 
because it was believed that corresponding contritions to 
the PNC amplitude are small and can be neglected. It 
was demonstrated in |12lll3l ] that radiative corrections to 
the PNC amplitude in cesium due to change of energy de- 
nominators and change of El transition amplitudes (see, 
eq. mSJ) are not so small when taken separately. But 
they cancel each other almost exactly when taken to- 
gether. 

It is interesting to note that very similar situation takes 
place for the Breit correction. As it was first demon- 
strated by Derevianko [|| and then confirmed by 
Dzuba et al 0, the effects of Breit interaction on the 
PNC amplitude in cesium due to change of energy de- 
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TABLE II: Relative Breit correction to weak matrix elements 
of Cs, Ba+, Tl, Fr, and Ra+. 



Atom/Ion s — p p — d 

Cs -0.58% -0.03% 

Ba+ -0.65% -0.29% 

Tl -1.17% -0.81% 

Fr -0.98% +1.60% 

Ra+ -1.03% +0.85% 



nominators and change of El transition amplitudes are 
not small, numerically very close, and have opposite 
signs. Corresponding numbers are -0.4% 

for the 

effect on El amplitudes and 0.3% [H and 0.4% for 
the effect on energies. There is almost exact numerical 
cancellation between two effects. Similar situation takes 
place for francium. 

However, as it is clear from the analysis of the s — d 
PNC amplitudes, numerical cancellation of the effect of 
Breit interaction on energies and El transition ampli- 
tudes is a fortunate feature of the s — s PNC amplitudes 
rather than a general situation. The same might be true 
for the radiative corrections to the PNC amplitudes. 

In Table [H] we present the relative values of the effect 
of Breit interaction on the s — p and p — d weak ma- 
trix elements between lowest valence states of five con- 
sidered atoms. These states give dominant contribution 
to the PNC amplitudes. Matrix elements are calculated 
in Brueckner approximation: 

W ab = (^ r \H w + SV w \^ r ). 

Note, however, that relative value of the Breit correction 
to the s — p weak matrix elements does not depend on 
the approximation. It is also the same for states with 
different principal quantum number n as far as the energy 
of this state is small compared to —Ze 2 /RN, where Rn 
is nuclear radius (see discussion above for more details). 
Breit correction to the PNC amplitude is approximately 
equal to the Breit correction to the weak matrix element 
only for the s — s PNC transitions of Cs and Fr. There is 
no obvious connection between the two for all other PNC 
amplitudes. This indicates a significant role of the Breit 
correction to the El amplitudes and the energies for the 
p — d and pi/ 2 — P3/2 PNC amplitudes (see eq. (tTC|l ). 

Another distinctive feature of the s — d PNC ampli- 
tudes which can be easily seen from Table H] is the signif- 
icant and non-trivial role of the correlations. One would 
expect that correlations either play no role as for the 
s — s PNC amplitudes or enhance the effect of Breit in- 
teraction. Indeed, correlations increase electron density 
on short distances which would lead to expectations that 
the role of the Breit interaction should also increase. In 
fact, the opposite happens. As it is evident from Table [i] 
correlations decrease the Breit correction to the s — d 
PNC amplitude. 

To understand this behavior let us start our analysis 
from single-electron matrix elements of weak interaction. 



Arguments in favor of the small role of correlations pre- 
sented above for the s — s PNC amplitudes were based 
on the fact that both Breit and weak interactions are lo- 
calized on short distances. While it is true for the "bare" 
operator of the weak interaction, it is not exactly true 
when many-body effects are also included. The effective 
operator of the weak interaction can be presented in the 
form 

H$ = Hw+SV w , (17) 

where first term is proportional to nuclear density and 
is zero everywhere outside the nucleus. Second term is 
due to change in self-consistent Hartree-Fock potential 
because of the effect of weak interaction. It describes the 
so-called core polarization effect and its radius is roughly 
equal to the radius of atomic core. For the s — p matrix 
elements core polarization is just small correction to the 
leading first term. Therefore all considerations based on 
assumption that weak interaction is localized on short 
distances are still approximately valid. 

For the p — d weak matrix elements first term in 11711 
gives zero for the point-like nucleus and all value for the 
matrix elements is given by the second, core polarization 
term. Even for finite size nucleus, contribution of the 
first term in l|17|) is several orders of magnitude smaller 
than of the second term. This means that p — d weak 
matrix elements don't come from short distances and all 
arguments based on that assumption are not valid for 
them. 

The EIpnc amplitude can be presented as a sum of 
three terms (see eq. (|1 5|) . Corresponding numbers for 
the s — d PNC amplitudes in Cs, Ba + , Fr and Ra + as 
well as for the 6pi/2 — 6P3/2 PNC amplitude in Tl are 
presented in Table ITTll The numbers in table are orga- 
nized in such a way that the term which is proportional 
to the s — p weak matrix clement goes first, second term 
is proportional to the p — d matrix element and the last 
term has the mix of both. Note that this last term is 
not included to the sum over single-electron states ex- 
pression (|16fl . It is small for alkali atoms. However, it is 
not negligible and we present it here to demonstrate that 
its inclusion is important for accurate calculations. 

The numbers in the table indicate that, as one would 
expect, term proportional to the s — p matrix element 
strongly dominates in the Elp^c amplitude. However, 
the term proportional to the p — d weak matrix element 
is not small at all and brings significant complications 
to the behavior of the EIpmc amplitudes. The case of 
particular interest is comparison of the s — d PNC am- 
plitudes of Cs and Ba + . As one can see from Tabic |T| 
the role of correlations in the Breit correction to these 
El pnc amplitudes is very much different in cases of Cs 
and Ba + . Correlations almost do not change Breit cor- 
rection to the s — d PNC amplitude in Ba + . On the other 
hand, their effect on the s — d PNC amplitude in Cs is sig- 
nificant. Moreover, it goes contrary to naive expectations 
as it was discussed above. The Breit correction is smaller 
when correlations are included while the effect of corre- 
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TABLE III: Effect of Breit interaction on different contributions to the s — d and p\n ~ P3/2 EIpnc amplitudes in Brueckner 
approximation (10 -11 ieasQw /N). 



Atom 


Transition 


Breit 


(V^l-Hfil + SV E l\5lpb) 


{tp a \H w + 5V w \X b ) 


{4>a\8V E \w\^b)- 


Sum 


or ion 


a — b 


(yes/no) 


~ s — p weak m.e. 


~ p — d weak m.e. 






Cs 


'-' u 'o/2 


no 


3.0045 


0.8314 


-0.1081 


3.7277 




yes 


2.9846 


0.8400 


-0.1081 


3.7165 






A 


-0.0199 


0.0086 


0.0000 


-0.0112 






A (%) 


-0.66 


1.0 


0.0 


-0.30 


Ba+ 


ou 3/2 u,b l/2 


no 


2.5776 


-0.2486 


-0.1465 


2.1826 




yes 


2.5583 


-0.2469 


-0.1462 


2.1653 






A 


-0.0193 


0.0017 


0.0003 


-0.0173 






A (%) 


-0.75 


0.07 


0.02 


-0.79 


Tl 


XJ l J Z/2 u^l/2 


no 


37.454 


-1.195 


-11.989 


24.270 




yes 


37.106 


-1.191 


-11.869 


24.046 






A 


-0.348 


0.004 


0.120 


-0.224 






A(%) 


-0.93 


0.03 


1.0 


-0.92 


Fr 


6d 3/2 - 7s 1/2 


no 


51.598 


6.466 


-1.627 


56.437 




yes 


51.074 


6.643 


-1.628 


56.090 






A 


-0.524 


0.177 


0.001 


-0.347 






A (%) 


-1.00 


2.74 


0.0 


-0.61 


Ra+ 


6^3/2 - 7Sl/ 2 


no 


46.767 


-1.964 


-2.294 


42.511 




yes 


46.230 


-1.973 


-2.287 


41.971 






A 


-0.537 


0.009 


0.007 


-0.540 






A (%) 


-1.15 


0.46 


0.31 


-1.27 



lations on the density of external electron on short dis- 
tances suggests that it should be rather larger. Another 
interesting thing is the difference between the results for 
Cs obtained in the S^ 2 ^ approximation (formula |Q ) 
and Brueckner approximation (formula Q15[0. These two 
expressions differ by higher orders in £. Therefore the 
difference in results indicates important role of high-order 
S terms. This is an interesting fact given that the term 
linear in £ strongly dominates in the correlation correc- 
tion to the PNC amplitude. Situation for the correlated 
Breit correction is opposite. The term linear in £ is small 
and higher-order terms dominate. 

For a more detailed analysis of the s — d PNC ampli- 
tudes in Cs and Ba + we present another set of data in 
Table HVI This data is similar to those in Table UTTl but 
in TDHF approximation. Data presented in Table ITTT1 is 
in Brueckner approximation which means that correla- 
tions are included to all orders in E. In contrast, data 
from Table llVl contains no correlations. Comparison of 
the data in Table IIIII and Table IIVI shows that correla- 
tion correction to terms proportional to the s—p weak 
matrix elements is small. This is similar to the s — s 
PNC amplitude and the same arguments are valid to ex- 
plain this. On the other hand, correlations significantly 
increase (by an order of magnitude) the contribution of 
the terms proportional to the p—d weak matrix elements. 
Correspondingly, Breit correction also increases. This is 
also well expected fact given that correlations increase 



density of external electron in the core. In summary, dif- 
ferent behavior of the s—p and p—d weak matrix elements 
with respect to correlations is due to different distances 
from which these matrix elements get their values. The 
s—p matrix elements are well localized at atomic nucleus, 
while p — d matrix elements localized on much larger dis- 
tances inside atomic core. 

Since correlations have different effect on Breit correc- 
tion to the 6S1/2 — 5^3/2 PNC amplitude in Cs depending 
on whether they are included in linear in E approxima- 
tion (formula dJ) of in all-order Bruekner approxima- 
tion (formula (|15fl ). it is important to check that all dif- 
ference is due to higher order in £ terms. To do this we 
extract linear in S contribution from Brueckner approx- 
imation and compare it with the result obtained using 
formula (|12l) . This is done with the use of rescaling pro- 
cedure similar to those discussed in the end of previous 
section for the Breit operator. However, this time we 
rescale correlation potential E. Namely, we replace E in 
eqs. l(Kl|) and ifT^I) by ££, where S is rescaling parameter. 
Calculating Breit correction with S < 1 and extrapolat- 
ing it to S = 1 leaves only linear in E contributions. The 
results are presented in Table W\ Extrapolation to S = 1 
is done using the formula 

A(S=l) = A(0) + a < S >: a <°>. 
Corresponding numbers are presented in the last column 
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TABLE IV: Effect of Breit interaction on different contributions to the 5d 3 / 2 — 6s PNC amplitude in Cs and Ba + in TDHF 
approximation (10 -11 ieasQw /N). 



Atom 
or ion 


Breit 
(yes/no) 


{i>Sd a/2 \HE1 + SVeI \5lp6s) 

~ s — p weak m.e. 


{ipsd 3/2 \Hw + 8Vw\X 6s ) 
~ p — d weak m.e. 


{•4>Sd 3/2 \$VElw\4>6a) ■ 


Sum 


Cs 


no 


2.9009 


0.2806 


-0.0637 


3.1177 




yes 


2.8836 


0.2814 


-0.0638 


3.1012 




A 


-0.0173 


0.0008 


0.0001 


-0.0165 




A(%) 


-0.60 


0.29 


0.0 


-0.53 


Ba+ 


no 


2.6079 


-0.3155 


-0.1245 


2.1678 




yes 


2.5886 


-0.3131 


-0.1242 


2.1513 




A 


-0.0193 


0.0024 


0.0003 


-0.0165 




A (%) 


-0.74 


0.76 


0.2 


-0.76 



TABLE V: Breit correction to the 6s — 5d$/2 PNC amplitude in Cs in Brueckner approximation as a function of rescaling 
parameter 5" (E -> SE in and JUJ) (10~ n ieasQw / 'N). 



S No Breit with Breit A A(ff) - A(0) A{S = 1) 

0.0 3.117696 3.101215 -0.016481 

0.1 3.160052 3.143462 -0.016590 -0.000109 -0.017571 

0.2 3.203363 3.186710 -0.016653 -0.000172 -0.017341 

0.4 3.295216 3.278633 -0.016583 -0.000102 -0.016736 

0.6 3.400763 3.384702 -0.016061 0.000420 -0.015781 

0.8 3.534769 3.520128 -0.014641 0.001840 -0.014181 

1.0 3.727607 3.716356 -0.011251 0.005230 -0.011251 



of Table One can see that the value of Breit correc- 
tion calculated at S = 0.1 and extrapolated to S = 1 
(-0.017571) practically coincide with the value -0.0176 
obtained using formula l|12[l (see Table P). For larger val- 
ues of S non-linear in E terms become important and 
rapidly change the value of Breit correction. 

There are few more important conclusions which can 
be drawn from the fact that Breit correction to a PNC 
amplitude is only sensitive to change of the electron den- 
sity. This change of density is proportional to the value 
of the correlation potential E. Therefore, it is important 
to check how inclusion of some extra terms into E affect 
the value of the Breit correction to the PNC amplitudes. 
For example, inclusion of the Breit interaction itself into 
calculation of E does not change its value very much. 
Corresponding change of removal energies of Cs is less 
than 0.02%, while change of the 6s — 6pi/ 2 energy interval 
is 0.04% only. Test calculations show that corresponding 
change to the PNC amplitude is also small, it is 0.04% 
for the 6s — 5g?3/2 PNC amplitude and even smaller for 
the 6s — 7s PNC amplitude. On the other hand, it is well 
known from a number of calculations (see, e.g. 0j) that 
inclusion into E of certain higher-order Coulomb terms 
leads to significant change in the value of E and electron 
density on short distances, brining theoretical values for 
energies, hyperfine structure, transition amplitudes, etc. 
to much better agreement with the experiment. There- 
fore, it would be inappropriate to assume a priori that 
higher-order Coulomb terms do not affect Breit correc- 
tion. This is to be checked numerically. We have per- 



formed such calculations and the results show that the 
effect of the higher-order Coulomb terms in E is indeed 
small. Its value is close to the effect of inclusion of Breit 
interaction into E. Moreover, the two effects tend to can- 
cel each other. One can say therefore that inclusion of 
Breit and higher-order Coulomb terms into E are equally 
important or unimportant depending of the needed ac- 
curacy. 

Let us now come back to comparison of the s — d PNC 
amplitude for Cs and Ba + for few more important obser- 
vations. If terms proportional to s — p and p — d weak 
matrix elements are compared separately, their behavior 
is very similar for both atoms. This is a natural conse- 
quence of the similar electron structure. However, the 
total s — d PNC amplitude is more sensitive to small dif- 
ference in the electron structure. In Ba + ion bd states 
are lower than the 6p states while in Cs 5d states are 
higher than the 6p states. This leads to larger correla- 
tion correction to the p — d part of the PNC amplitude 
in Ba + compared to those of Cs. This also leads to the 
opposite sign of these terms in Ba + and Cs. As a result, 
the effect of correlations on final PNC amplitude is very 
different. In general, the effect of correlations on s — d 
and pi/2 — Ps/2 amplitudes is very hard to predict due 
to different behavior of terms proportional to s — p and 
p — d weak matrix elements. 

Having accurate results for Breit correction to the PNC 
amplitudes for many atoms we can now check how they 
rescale with the nuclear charge Z. This would allow us 
to obtain Breit correction for other atoms without doing 
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TABLE VI: Breit corrections to the PNC amplitudes in Cs, 
Ba + , Tl, Fr and Ra + (10 -11 ieasQw / N), comparison with 
other calculations. 



Atom 


Transition 




Source 


Cs 


6s — 7s 


-0.0056 


this work 






-0.0054 


Derevianko ^4J 






-0.0055 


Dzuba et al ^ 






-0.004 


Kozlov et al \8j 






-0.0045 


Shabaev et al [21J 




6s - 5d 3/2 


-0.0112 


this work 


Ba+ 


6s - 5d 3/2 


-0.0172 


this work 


Tl 


6P1/2 - 6p 3/2 


-0.224 


this work 


Fr 


7s - 8s 


-0.15 


this work 






-0.18 


Derevianko 






-0.15 


Safronova 








and Johnson [20J 






-0.14 


Shabaev et al ^lj 


Ra+ 


7s - 6d 3/2 


-0.541 


this work 



sophisticated calculations. It is easy to see that Breit 
correction to similar amplitudes in atoms with similar 
electron structure is just proportional to Z. For exam- 
ple, AElpjyc(Fr) = (87/55) A J Bl PA r C .(Cs). It is inter- 
esting to note that rescaling from the 6s — 7s amplitude 
in Cs to the 6p\/2 — &P3/2 amplitude in Tl gives result 
which is very close to the correct answer: AiJlpjvc(Tl) = 
(81/55) AEl PNC (Cs) = -0.91% while accurate calcula- 
tions give —0.92%. However, this is just fortunate coin- 
cidence. Separate analysis of terms proportional to the 
s — p and p — d weak matrix elements show that no reli- 
able rescaling between Cs and Tl is possible. The same 
is true for any other pair of the PNC amplitudes which 
correspond to different electron structures. 

We present our final results in Table IvTl together with 
the results of other authors. The results for the 6s — 
7s amplitude in Cs and 7s — 8s amplitude in Fr are in 
good agreement with other calculations. However, we 
believe that our results are the most accurate due to non- 
perturbative treatment of the Breit interaction which is 
complete up to second order of the MBPT. The results 
for other atoms as well as for the s — d PNC amplitudes 
for Cs and Fr are presented for the first time. 

We now can use the result of this work for the Breit 
interaction contribution and our previous calculations of 
other contributions to find the PNC amplitude in 133 Cs, 
extract value of the weak charge and compare it with 
the standard model. Our many-body calculations of the 
PNC amplitude produced by the electron-nucleus weak 
interaction gave the following result [Tjj (without the 
Breit contribution and QED radiative corrections) 
E PNC = -0.9060(1 ± 0.5%) x lQ- n iea B (-Qw/N) . 
We take the Breit contribution -0.0056 from this work 
and the QED radiative corrections -0.0029 (-0.32 % of 
the PNC amplitude) from Ref. 12]. This gives 

E PNC = -0.8975(1 ±0.5%) x W- n iea B (-Qw /N) . 

_ (18) 
In the experiment [2j the ratio Ep^c/P, where (3 is the 



vector transition polarizability, actually has been mea- 
sured. To extract Q w we used (3 = 26.99(5)a| obtained 
in 0| as the statistical average of the two most accu- 
rate values of (3. From the measurements of the PNC 
amplitude Q we obtain 

Qw = -72.69(29) cxp (36)thcor . (19) 

The difference with the standard model value Qyy 1 = 
-73.19(3) |3 is 

Qw ~ = 0.50(46) , (20) 

adding the errors in quadrature. Thus, the difference 
between the central points is one standard deviation. 

For 205 T1 we use PNC amplitude from Ref. 16], 
E PNC = -(27.0 ± 0.8) x lQ- n iea B (-Qw/N), 
the QED radiative correction -0.57 % (see detailes in |]|), 
the neutron skin correction -0.3 % and the Breit cor- 
rection -0.92 % from this work. The result is 

E PNC = -(26.5 ± 0.8) x W- n iea B (-Qw/N) . (21) 

Using Ml amplitude 1.693 a.u. from |23j and measure- 
ment Im(E PNC /Ml) = -(14.68 ±0.17) x 10" 8 a.u. 
we obtain for 205 Tl 

Q W = -116.2(1.3)exp(3.5) the or . (22) 

The standard model value is -116.81(4) The differ- 
ence with the standard model for Tl is 

Qw -Q S w= 0.6(3.7) , (23) 

IV. CONCLUSION 

We have developed a method of non-perturbative 
treatment of the Breit interaction in accurate calcula- 
tions for many-electron atoms. The method has been 
applied to calculation of the Breit correction to the PNC 
amplitudes of Cs, Ba + , Tl, Fr and Ra + . The results for 
the 6s — 7s amplitude in Cs and 7s — 8s amplitude in Fr 
are in good agreement with other calculations. The re- 
sults for other atoms and for the s — d PNC amplitudes for 
Cs and Fr are presented for the first time. The following 
features of the Breit correction to the PNC amplitudes 
have been revealed: 

• The effect of correlations on the s — s PNC ampli- 
tudes is small. 

• In contrast, correlations are important for those 
PNC amplitudes which depend on the p — d weak 
matrix elements between valence states. 

• For accurate treatment of correlations inclusion of 
linear in the correlation potential E terms is not 
enough. The main effect is due to the change in 
density of external electron inside the core which is 
magnified in higher orders in E. It can be included 
by the use of Brucckncr orbitals. 
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• Inclusion of the Brcit interaction as well as higher- 
orders Coulomb terms into correlation potential £ 
have only small effect on PNC amplitudes and can 
be neglected in most cases. 

• Numerical cancellation of the effects of Breit in- 
teraction on El transition amplitudes and energy 
denominators in the s — s PNC amplitudes for Cs 
and Fr is rather fortunate. No such cancellation 
takes place for other atoms and for the s — d ampli- 
tudes of Cs and Fr. The same is probably true for 
the radiative corrections to the PNC amplitudes. 



• Rescaling of the Brcit correction to the PNC am- 
plitudes using ratio of the nuclear charges Z can 
only be done to similar amplitudes in atoms with 
similar electron structure. 



Combining Breit corrections to PNC amplitudes in ce- 
sium and thallium with previous accurate calculations 
and measurements confirms good agreement of the weak 
nuclear charges of these atoms with the prediction of the 
standard model. 
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